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Abstract 

O 

The renormalization of the weak charge-changing hadronic current as a func- 
tion of the reaction energy release is studied at the nucleonic level. We have 
calculated the average quenching factors for each type of current (vector, axial 
vector and induced pseudoscalar). The obtained quenching in the axial vec- 

X" 

' tor part is, at zero momentum transfer, 19% for the sd shell and 23% in the 

fp shell. We have extended the calculations also to heavier systems such as 
56 Ni and 100 Sn, where we obtain stronger quenchings, 44% and 59%, respec- 
tively. Gamow-Teller type transitions are discussed, along with the higher 
order matrix elements. The quenching factors are constant up to roughly 60 
MeV momentum transfer. Therefore the use of energy-independent quench- 
ing factors in beta decay is justified. We also found that going beyond the 
zeroth and first order operators (in inverse nucleon mass) does not give any 
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substantial contribution. The extracted renormalization to the ratio Cp/Ca 
at q = 100 MeV is -3.5%, -7.1%, -28.6%, and +8.7% for mass 16, 40, 56, 
and 100, respectively. 
PACS number(s): 21.30.Fe, 23.40.-s, 21.60.Cs, 21.10.Pc 

I. INTRODUCTION 

The phenomenological structure of the weak hadronic current between the proton and 
neutron states is well determined by its properties under the Lorentz transformation. The 
additional constraints come from the requirement of time reversal symmetry as well as 
from the invariance under the G-parity transformation (combined charge conjugation and 
isospin rotation). The resulting interaction Hamiltonian consists of vector (V), axial vector 
(A), induced weak magnetism (M), and induced pseudoscalar (P) terms together with the 
associated form factors C a , a = V, A, M, or P. These form factors are called as coupling 
constants at zero momentum transfer. The present experimental knowledge does not exclude 
the presence of the scalar and tensor interactions. However, their contribution is expected 
to be small due to weak couplings Q. 

The values of vector, axial vector, and weak magnetism couplings are well established 
by beta-decay experiments as well as by the conserved vector current hypothesis (CVC), 
introduced already in the late 50 's @. The magnitude of the pseudoscalar coupling is more 
uncertain, although the partially-conserved axial current hypothesis (PCAC) provides an 
estimate along with muon capture experiments in hydrogen [|]||. The value of Cp in the 
nuclear medium is not precisely established. 

In nuclear beta decay, with an energy release up to some 20 MeV, only the vector (Fermi) 
and the axial vector (Gamow-Teller) terms are usually important. The induced pseudoscalar 
and weak magnetism parts are essentially inactive, since their contributions are proportional 
to q/M, where q is the energy release and M is the nucleon mass (in units where H = c = 1). 
There are, however, weak nuclear processes, like muon capture, where the energy release is 
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much higher (in muon capture, typically q ~ m M pa 100 MeV). 

Summed theoretical beta decay strengths are systematically larger than the experimental 
ones. This so-called quenching of the (allowed) decay strength is usually explained in terms 
of core polarization (degrees of freedom which are left out from the model space) and non- 
nucleonic degrees of freedom like isobars and meson exchange currents ||. Many authors 
(e.g., |6|-f|) have established the quenching factors for the Gamow-Teller decays and closely 
related magnetic dipole (Ml) transitions. 



In [|lQ| , [liri we have self-consistently constructed effective operators for the weak hadronic 
current between proton and neutron states. These operators, as explained in Sec. |T|, take 
into account the above-mentioned core polarization effects, which are expected to be the 



largest correction to the bare matrix elements [fL2| . In the present work our aim is to 



calculate self-consistently the quenching for all types of operators (V, A, M, and P) for 
energies up to the muon capture range. In addition to the traditional shell-model regimes, 
sd and fp shells, we have extended our calculations to |gNi and ^[jSn as closed-shell cores. 

Typically one expects about 20% quenching of the axial vector part in the sd shell, i.e., 
the calculated Gamow-Teller matrix elements (a) are to be multiplied by a factor ~ 0.8 || 
when both core polarization and non-nucleonic degrees are accounted forQ Effects of the 
same magnitude are expected in the fp shell ||. The authors in |J reach the conclusion that 
the mass dependence of the quenching has saturated already in the fp shell. However, the 
major shell closures which separate the spin-orbit partners in 56 Ni (Of 7/2 and 0/ 5 / 2 , 0g g / 2 and 
0^7/2) and 100 Sn (0g 9 / 2 and 0g 7 / 2 , 0h u / 2 and 0h 9 / 2 ) model spaces introduce large first-order 
corrections to the operators. Thus, the situation is not analogous to the one seen in light 



nuclei with closed LS shells. For a recent work in the Sn mass region, see e.g., Ref. p^] . We 
remind the reader that the quenching is always related to the choice of the model space. 
The nuclear muon capture process can be used to extract the ratio Cp/Ca- Unfor- 



1 The Gamow-Teller strengths -B(GT) oc (a) 2 are then multiplied by (0.8) 2 ~ 0.6. 
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tunately, the results for partial capture rates are very sensitive to applied nuclear model, 



especially to the model space and to the residual two-body interaction (see, e.g., |T3[ and 
references therein). The total rates offer perhaps a more reliable source of information, in- 
dicating no or only small quenching for the ratio Cp/Ca ||15| . It is of interest to see whether 
this quenching can be explained in terms of effective charges for the axial vector and pseu- 
doscalar operators, that is, without the complications coming from the nuclear structure 
calculations. 

In addition to the zeroth-order Fermi and Gamow-Teller type operators, our set includes 
the first-order terms in the transition amplitude (first order in q/M as well as velocity- 
dependent terms). We shall also examine the importance of the second-order terms. We 
stress that the results obtained in this work can be applied quite generally. We have used 
the muon as an initial bound state lepton, but the results apply to electron capture as well 
and therefore to beta decay in general (in our calculations, the muon is nothing but a heavy 
electron!). 

In this work, after a short review of the formalism of the semileptonic weak processes 
in Sec. |TJ and effective operators in Sec. fTJ, we concentrate on the results in Sec. |TV|. We 
consider four cases, with 16 0, 40 Ca, 56 Ni, and 100 Sn as closed-shell cores. 

II. INVARIANT AMPLITUDE AND SINGLE-PARTICLE OPERATORS 

After the standard nonrelativistic reduction, the semileptonic charge- changing weak pro- 
cess 

X b +p^ux + n, (1) 

where A& is a bound (anti)lepton in an atomic 15* orbit and V\ is the corresponding 
(anti)neutrino, is described by the amplitude 

M 2 = Y j \Mv{^u) + M a {k,u) + M p {k,u)\ 2 . (2) 
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We take to be a muon, with a mass m M = 105.658 MeV. The form of the effective weak 
hadronic current used for Eq. (0) is the most general one consistent with the expected G- 
parity symmetry and time reversal symmetry. The functions M a (n,u) include the form 
factors, transition operators, and angular momentum couplings. Explicitly, the vector part 
is given by 

M v (K,u)/C v (q 2 ) = [01u]S Ou (k)5i u - ^[lhp]S' lu (- K ) 



+ ^/f^t 1 +fh~ Mn) {y/f+lWillul, 11 + l)[ll+ lu+] (3) 
+ \fiw(llul 7 ll- l)[ll- Iu-]\s' 1u (-k). 

We consider the V and M terms together, as suggested by the CVC: Cm = (fMp— }J, n )Cy/2M, 
where fi p and fi n are the anomalous magnetic moments of the proton and neutron in nuclear 
magnetons. The Fermi- type term is the first term on the right hand side of Eq. (^). The 
axial vector part is given by 

M A (K,u)/C A (q 2 ) = -[11u]S 1u {k) + jj[0lup]5 lu S' 0u {-K) - M P (n,u)/C P (q 2 ), (4) 

including the pseudoscalar part 



M P (K } u)/C P (q 2 ) = - 7r ^ 77 { t/irWlli + 1«+] + 1u-]\5 Iu S' 0u (-k). (5) 



l + l 

2VSM \ V 21 + 1 



The Gamow-Teller-type term is the first term on the right hand side of Eq. @ . 
In Eqs. k labels the quantum numbers of the emitted neutrino v\, 

k>0: j = l-- 1 = k (6) 
K < : j = I + i I = -k - 1, (7) 

where I and j are the orbital and total angular momentum quantum numbers of v\. The 
quantity I is given by I — Sign(«), W are the usual Racah coefficients, and 
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Sku(K) = J2(2j + l)W(l/21jl,l/2u), k = l 



S L(- K ) = Sign(/t)S fcu (-/t). 

The most important ingredients, the transition operators, are embedded in the reduced 
matrix elements [/cwit], [kwu±], and [kwup]. The quantum numbers labelling the matrix 
elements are k = S\ + s u , so that k = |k| = or 1, and u = |u| = |J/ + Jj| is the tensorial 
rank of the transition operator. The symbol w is the rank of the spherical harmonics and 
is therefore related to the parity change. It is given by w = I for [kwu] and [Irwui] type 
matrix elements, and w = l + l or w = l — 1 for [kwup] type matrix element (k and w 
must be able to couple to u). The symbol p labels the momentum dependent operators. 
The matrix elements with the corresponding single-particle operators are listed in Table |. 
These operators are further multiplied by the radial wave function of the initial state lepton 
16| . We have taken into account the large component 



G M (r) = 2(aZm' fl ) 3 / 2 e- aZ < r , (9) 

where a ~ 1/137 is the fine structure constant and m' is the reduced muon mass. 

The amplitude (|2]) can be used for the calculations of muon (or electron) capture rates 
]T^,|TB|. As mentioned in the previous Section, our aim is to calculate the effective charges 
(effective form factors) for the vector, axial vector, and pseudoscalar parts of the amplitude, 
so as to help to understand the differences between calculated and experimental rates and 
other observables. 

For the actual calculations, we divide the reduced nuclear matrix elements into single- 
particle and many-body parts, 

M a (K, u) = 5>| \m a (K, u) | \p) [JM ^ p] y i \ (10) 

pn V^J + t 

where n = (n n , l n ,j n ) and p = (n p , l p ,j p ) label the single-particle states. The doubly-barred 
matrix elements are reduced in the angular momentum space. These reduced single-particle 
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matrix elements, which we calculate in the harmonic oscillator basis with the single-particle 
operators m a (K,u), are constructed as described in the following Section. Further, we have 
defined 



a jm = (-iy +m a 



'3,-m- 



(11) 



As an example, from Eq. (01), the single-particle matrix element for the axial vector part 



is 



(n\\m A (K, u)\\p) = -C ' A (q 2 ){n\\Ou u G At (r)||p)S , i n (/t) + 



C A (q 2 
M 



n\\O oI G^r) | \p)S Iu S' 0u (-k) 



+ C A (q 



2 n q 



l + l 



2V3M V 2/ + 1 



[n\\°iT+iu+Gi*(r)\\p) 



(12) 



+ 



I 



21 + 1 



{n\\°ii-iu-G^(r)\\p) } $iuS'ou(- K )- 



A closer look into the Gamow-Teller type matrix element gives 



CU(g 2 )(n||CW^(r)|b)£i„(K) = -C A ( q y^(l nJn \\Y Uu (r,a)\\l p j p ) 



x / r 2 R nnln (r)3Mr)R nplp G,(r)drp(2 3 + 1) W{\/2\jl, 1/2 u). (13) 



Here R n i{f) are the radial single-particle wave functions. The reduced matrix element of the 
vector spherical harmonics is given by 



(lnjn\\Yu u (i,a)\\l p j p ) 



167T 



;(-l) lp+jp+jn+l+1 luj p j r 



( . . \ 

Jp Jn U 



1/2 -1/2 



(14) 



2u(u+l) 



i i\ 



\ 



u 

1 -1 



+ (-1) 



l n + l/2+j n 



u 





where x = \j2x + 1. The expressions for the vector and pseudoscalar parts are obtained in 
a similar way using Eqs. @ and @ and Table [I], where the operators O of Eq. (|1^) are 
given. For more details, see |L6| . 

The many-body part, which includes the one-body transition density (OBTD), is given 
by the adopted nuclear model. The effective operators introduced in Sec. [HI] do not affect 



the OBTD, which are taken to be given numbers. We do not calculate them here (see, e.g., 
fl4| for examples). In what follows we shall consider only the single-particle part of Eq. 
(|T0|). This is analogous to, e.g., Ml-transitions, where the effective charges are calculated 
without referring to the many-body part. 



III. PERTURB ATIVE METHODS AND EFFECTIVE OPERATORS TO SECOND 

ORDER 

In order to derive a microscopic approach to the effective operator within the framework 
of perturbation theory, we need to introduce various notations and definitions pertinent to 
the methods exposed. In this section we briefly review how to calculate an effective one- 
body operator within the framework of degenerate Rayleigh-Schrodinger (RS) perturbation 
theory |TT| , [T8|| , see also Refs. P,p~2f for a detailed discussion on various effective operator 



diagrams. 

It is common practice in perturbation theory to reduce the infinitely many degrees of 
freedom of the Hilbert space to those represented by a physically motivated subspace, the 
model space. In such truncations of the Hilbert space, the notions of a projection operator 
P onto the model space and its complement Q are introduced. The projection operators 
defining the model and excluded spaces are defined by 

p = j2mm, (15) 
i=i 

and 

Q= J2 (16) 

i=D+l 

with D being the dimension of the model space, and PQ = 0, P 2 = P, Q 2 = Q, and 
P + Q — I. The wave functions are eigenfunctions of the unperturbed hamiltonian 
Hq = T+U, where T is the kinetic energy and U an appropriately chosen one-body potential, 
that of the harmonic oscillator (h.o.) in this calculation. The full Hamiltonian is then 



S 



rewritten as H = H + Hi with Hi = V — U, V being e.g., the nucleon-nucleon (NN) 
interaction or the G-matrix to be discussed below. The eigenvalues and eigenfunctions of 
the full Hamiltonian are denoted by \^ a ) and E a , i.e., 

H\* a ) = E a \* a ). (17) 

Rather than solving the full Schrodinger equation above, one defines an effective Hamiltonian 
acting within the model space such that 

PH cS P |* a ) = E a P |* a ) = E a |$ a ) (18) 

where = P |\I/ a ) is the projection of the full wave function onto the model space, the 
model space wave function. In RS perturbation theory, the effective interaction H e g can be 
written out order by order in the interaction Hi as 

PH cS P = PHiP + PHi—HiP + PHi—Hi—HiP + .... (19) 

e e e 

Here we have defined e = uo — H , where oo is the so-called starting energy, defined as the 
unperturbed energy of the interacting particles. Similarly, the exact wave function \^f a ) can 
now be written in terms of the model space wave function as 

|*„> = |*«> + -Hi |<U + ^-H^-Hi |$ a > + . . . (20) 
e e e 

In studies of nuclear transitions such as beta decay, the quantity of interest is the tran- 
sition matrix element between an initial state and a final state of an operator O 
defined as 

= (*w . (21) 

Since we perform our calculation in a reduced space, the exact wave functions \^f t i) are 
not known, only their projections onto the model space. We are then confronted with the 
problem of how to evaluate Ofi when only the model space wave functions are known. In 
treating this problem, it is usual to introduce an effective operator O^, defined by requiring 
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O fi = (^ f \O eS \^). (22) 

Observe that O e s is different from the original operator O^. The standard empirical proce- 
dure is then to introduce some adjustable parameters in O^. 

The perturbative expansion for the effective operator can then be written in a much 
similar way as Eqs. flT9|) and (]2I1|), i.e., 

(23) 

In Fig. [I] we list all diagrams (except folded diagrams) to second order in the interaction 
evaluated in this work. We do not include Hartree-Fock insertions. For pure Gamow-Teller 
or Fermi like operators, see e.g., the review article by Towner ||, such diagrams are exactly 
zero. Another feature of e.g., the Gamow-Teller-type operators is that for several diagrams 
involving particle-hole contributions, these diagrams are exactly zero unless the particle-hole 
orbits are spin-orbit partners. This means that for LS-closed shell nuclei like 16 O and 40 Ca 
diagrams like (I)-(VIII) or (XIII)-(XX) are all zero. However, this picture changes when we 
move to closed-shell nuclei like 56 Ni and 100 Sn. For Ni the last proton and neutron holes 
are in the Of 7/2 single-particle orbit. This means that the Of 7/2 hole and the O/5/2 particle 
states in 56 Ni yield non-vanishing contributions to the Gamow-Teller-type operator from the 
above-mentioned diagrams. Similarly, in Sn these contributions are represented by the spin- 
orbit partners in the the 0*79/2 hole and the 0*77/2 particle states. These spin-orbit partners 
yield then lhu intermediates states. Similarly, we have also spin-orbit partners for particles 
states outside the model space. These are 0*79/2 and 0*77/2 for 56 Ni and Ohn/2 and Ohg/2 for 
100 Sn. These lhu intermediates states are then responsible for the different quenching of 
the effective operators in the mass regions of 16 O- 40 Ca and 56 Ni- 100 Sn, respectively. 

We end this section with a discussion of how to construct a G-matrix. The G-matrix 
enters in turn our perturbative expansion for the effective operator. As is well know in 
nuclear physics, the NN potential exhibits a repulsive core, which renders any perturbative 
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treatment prohibitive. However, one possible way of overcoming this deficiency is to intro- 
duce the reaction matrix G, which accounts for short range correlations. The G-matrix is 
defined through 

G = V + V %—G. (24) 

uj - QTQ 

Here, u is the energy of the interacting nucleons in a medium and V is the free NN potential. 
We have assumed that the energy of the intermediate states can be replaced by the free 
kinetic spectrum T, since these states are predominantly of high excitation energy. 

In this work we solve Eq. fl24|) for finite nuclei by employing a formally exact technique 
for handling Q, originally presented by Tsai and Kuo JjlJ and discussed in Ref. [|(J. Tsai 
and Kuo employed the matrix identity 

Q QlQ Q= l-J P php P l' (25) 



with A = to — T — V, to rewrite Eq. (|24|) as 

G = G F + AG, (26) 

where Gp is the free G- matrix defined as 

G F = V + V—^G F . (27) 

The term AG is a correction term defined entirely within the model space P and given by 

AG = -V-P P^V. (28) 

A PA~ l P A v ; 

Employing the definition for the free G-matrix of Eq. (|27|) , one can rewrite the latter equation 



as 

AG = -G F -P—- l — —P-G F , (29) 

e P^e- 1 + e~ 1 G F e- 1 )P e ' V ; 

with e = to — T. 

We see then that the G-matrix for finite nuclei is expressed as the sum of two terms; 
the first term is the free G-matrix with no Pauli corrections included, while the second term 
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accounts for medium modifications due to the Pauli principle. The second term can easily 
be obtained by some simple matrix operations involving the model-space matrix P only. 

Finally, in order to calculate the G-matrix for the various mass regions, we need to define 
the relevant model spaces used to define the P and Q operators in the equation for G. The 
oscillator energies TlVL will be derived from TiVL = 45y4 -1 / 3 — 25A~ 2 / 3 , A being the mass 
number. This yields hQ = 13.9, TiVL = 11.0, HQ = 10.05, and HQ = 8.5 MeV for A = 16, 
A = 40, A = 56, and A = 100, respectively. We choose the model spaces which are believed, 
from both experiment and theoretical calculations, to be relevant as a first approximation 
for calculations of effective interactions and operators in the mass areas from A = 16 to 
A = 100. These are the 0d 5 / 2 , 0d 3 / 2 , and ls x / 2 orbits for A = 16, the lp3/ 2 , lpi/2, O/7/2, 
and O/5/2 orbits for nuclei in the mass region of A = 40, the 1^3/2, lpi/2, 0/ 5 / 2 , and 0g 9 / 2 
orbits for nuclei in the mass region of A = 56 and the 0hn/ 2 , 0g 7 / 2 , ld 5 / 2 , ld 3 / 2 , and 2si/ 2 
orbits for A = 100. For these systems, the closed-shell cores ( 16 0, 40 Ca, 56 Ni, and 100 Sn) 
have equal numbers of protons and neutrons, and the model spaces are the same for both 
protons and neutrons. 

The definition of the Pauli operator for the G-matrix can be found in Refs. P0"|,|2T 



where the so-called double-partitioned scheme has been used. A detailed discussion on the 



computation of the G-matrix can be found in Ref. This definition means that also the 
shell above that which defines the model space of the effective interaction, is included in the 
evaluation of the G-matrix. For the IsOd-shell, this means that we also include the lpOf shell 
in the definition of the P-operator for the G-matrix. As a consequence, we have to include in 
our perturbation expansion ladder type of diagrams where the allowed intermediate states 
are those of the lpO /-shell. With this prescription, we have to evaluate diagrams X, XIII, 
XIV, and XXIII-XXVIII in Fig. 0. 

In our actual calculation of the various effective operators, we truncate the sum over 
intermediate states at excitations of 4 — 8hu> in oscillator energy. This truncation yields and 
error of ~ 1% in our evaluation of the effective operator. The nucleon-nucleon interaction 



employed in this work is the CD-Bonn interaction of Machleidt et al. |22 
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IV. RESULTS 



We consider the transition + — > 1 + , which is of rank u = 1. In addition, k = —1 or 
k = 2, corresponding to j = | and / = or j = | and / = 2, respectively. The nuclear matrix 
elements, allowed by these quantum numbers together with parity conservation, are [101], 
[121], [101—], [121+], [Ollp], and [lllp]. We remind the reader that the matrix element [101] 
is closely related to the Gamow-Teller matrix element of the nuclear beta decay, only the 
radial dependence is more complicated due to the possibility of a larger energy release (see 
Table |). Our set includes matrix elements which are classified as forbidden in the nuclear 
beta decay |L6 . 



In order to discuss the average quenching, we define a factor p c 

E P „ \(n\\m™ n (K,u) 



Pa{q) - pa -El\(n\\m^( K ,u)\\ P W m 
where a = V, A, or P, and "ren" and "bare" refer to renormalized and bare single-particle 

matrix elements, respectively. The sums run over all the single-particle states included 

in the model space. We use absolute values in the sums in order to avoid cases where two 

single-particle matrix elements, with same magnitudes and opposite signs, cancel each other. 

This kind of cancellations do not easily happen in nuclear structure calculations, since the 

involved OBTD have different magnitudes (and signs). 

We start the discussion with sd shell nuclei. The effective operators are calculated with 

X gO as a closed-shell core. The model space is the full IsOd shell. From Fig. ^|we see that the 

quenching for V, A, and P terms remains essentially constant for the whole energy range 

considered. In particular, at the beta decay energies (below 20 MeV) all k = — 1 terms are 

constant. Moreover, we have p\(0) = 0.81 which is somewhat larger value than the empirical 

"universal quenching factor" of Ref. |8[], fitted to a large body of beta decay data in the sd 

shell (we remind the reader that our effective operators do not include the subnucleonic 

degrees of freedom). We also get a clear renormalization in the vector and pseudoscalar 

parts of the current. Note that CVC does not apply here, since we are not looking at Fermi 

transitions, the vector-type contribution comes from the higher-order terms. 
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Strictly speaking, our quenching factors are applicable only for one-particle systems, e.g., 
for 17 0. In practice, these factors are used for the whole model space, and the (weak) mass 
dependence is simply left out. In fact, we find that p\(0) = 0.81 for ^fSi as well (at the one- 
particle level). At the end of the model space, effective operators derived for the hole state 
can also be used. Then, often a larger renormalization compared to the particle operator is 
seen [|], reflecting the mass dependence. 

For 2oCa as a core, Fig. |3|, the overall features are very similar to 16 0. We get slightly 
more quenching, p\(0) = 0.77. This is in line with earlier studies, where only a slightly 
larger quenching for the fp shell is introduced. Our model space includes all the single- 
particle orbits of the Oflp shell. In || the authors reach the conclusion that already in the 
fp shell the quenching factor has reached the large- A limit. This is, indeed, confirmed by 



our results. We also note the good agreement with the results of Ref. ||23|| , Table 1. 

Qualitatively, the saturation of the quenching comes from the similar choice of the model 
space (complete Ohu space) in 16 and 40 Ca. In both cases, the first-order diagrams give zero 
contribution. Therefore the second-order contribution, with lhu intermediate excitations, 
is the most important one, and very similar behaviour is expected and, indeed, seen. The 
differences can be attributed e.g., to different oscillator parameters and differences in the 
single-particle orbit structures (in 40 Ca more ITiu excitations are available). Therefore, 
whenever the model space includes the whole Ohu oscillator shell, a similar quenching is to 
be expected. The small variations depend on the mass of the closed-shell core. 

The model space in which the effective operators are calculated for ggNi does not have 
a closed LS core. Now, the first-order transitions between / 7 / 2 , ^5/2, #9/2, and #7/2 single- 
particle orbits become possible (the situation is analogous to the Ml operator, which is 
diagonal in orbital angular momentum and spin). This is clearly reflected in the values 
shown in Fig. We now have p\{0) = 0.56. We also remind the reader that the Ikeda sum 
rule for Gamow-Teller beta decays is not fulfilled in this space. 

At mass A = 100 our space includes the single-particle orbits 0^7/2, ld 5 / 2 , ld 3 / 2 , 2si/ 2 , 
and 0/in/2 above the N = Z = 50 (^[jSn) major shell closure. The spin-orbit partners of 
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0(?7/2 and Ohn/2 orbits, 0g 9/ / 2 and 0h 9 / 2 , are missing from our model space. As in 56 Ni, this 
is seen from Fig. ^|: All terms are quenched by a factor which is clearly larger than in lighter 
nuclei with closed LS shells. In the axial vector part we have now p\(0) = 0.41, which again 
stays nearly constant up to tens of MeV's, well beyond the beta decay energy range (see 



also the discussion in Ref. JIB] ). The diagonal transitions g 7 / 2 — > #7/2 and /in/2 — > ^11/2 are 
responsible for the kink between q = 80 and 100 MeV. 

At the single-particle transition level, the spin-flip matrix elements (e.g., fj/ 2 — > $5/2) 
are evidently more quenched than the diagonal ones. This feature does not depend on the 
operator or mass. We further note that the Gamow-Teller type matrix element [101], being 
by far the dominant axial term, follows very closely the trend of the axial vector part through 
the whole mass range. 

The k = 2 axial vector terms are shown in Fig. |6|. Only the axial vector part is shown, 
since V and P terms are identical to the 1. For A = 16, 40, and 56 a strong 

quenching is obtained at high momentum transfers. This is mainly caused by the spherical 
Bessel function j w (qr). Oscillations in pa{q) are a sign of the interference between j w (qr) and 
the radial single-particle wave function. Clearly, a coherent extraction of the renormalization 
is not as feasible as in the k = — 1 case. 



In the report by Ciechanowicz et al. [^J] the meson exchange contribution to the muon 
capture matrix elements was found to be very small, at least in capture by 28 Si. This is, 



however, in contrast with the results for A = 12 nuclei of Ref. ||25|| . The quenching of 
the spin matrix element, essentially the Gamow-Teller matrix element, is expected to be 
dominated by the core polarization correction |J. However, our results leave some room 
for the subnucleonic corrections. For example, in the sd shell, about 50% of the observed 
quenching |J comes from the A isobars, meson exchange currents and more complicated 
many-body terms. The situation is similar also in the beginning of the fp shell ( 40 Ca). 

The second order terms in the invariant amplitude describing the process (|1|) are propor- 
tional to M~ 2 . Therefore one would expect their contribution to be at most a few percent of 
the first order terms [see Eqs. (§), (0), and @]. Indeed this is the case. For example in 16 0, 
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the r-dependent single-particle matrix elements which, according to Barabanov ||26|| , should 
be the dominating second order terms are about order of magnitude smaller than the first 
order terms. Only seldomly are the magnitudes roughly equal (one case in 16 0). Thus, when 
the first-order terms are scaled by M~ l and second-order terms by M~ 2 , the magnitudes 
behave roughly as 50 : 1, respectively. Their contribution is, in this context, neglible. The 
same conclusion is reached in ||26|| , where detailed expressions for the matrix elements are 
given. 

The quenching factors pa and pp can be used to estimate the ratio Cp/Ca- We take 
the data from Figs. |2|, |3], f|, and || at q = 100 MeV, which corresponds to the muon capture 
region. Then we have —3.5%, —7.1%, —28.6%, and +8.7% changes in Cp/Ca for masses 
16, 40, 56, and 100, respectively. If the bare value is taken from PCAC, Cp/Ca ~ 8.4, we 
have Cp/Ca ~ 8-1? 7.8, 6.0, and 9.1 for masses 16, 40, 56, and 100. This yields an average 
~ 7.8. Although not directly comparable to our results, it is interesting to note the results 
of Kolbe, Langanke, and Vogel |15[]. They used continuum random phase approximation 
to calculate the part of the capture rate which goes above the particle emission treshold. 
Their results show a reasonable agreement with data when the bare couplings are used. This 
nicely demonstrates the fact that if the model space dimension is increased, the couplings 
should asymptotically reach the bare values. 



V. SUMMARY 

We have constructed the effective transition operators corresponding to the general form 
of the weak hadronic current between the proton and neutron states. The effects of the 
renormalization are investigated as a function of the transition g-value, and an average over 
the single-particle transitions is taken separately for vector, axial vector and pseudoscalar 
terms. We have considered only nucleonic degrees of freedom. In addition to the operators 
present in the allowed beta decay, we have considered the higher-order corrections to the 
transition amplitudes. 
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In the sd and fp shells, we get 19% and 23% quenchings in the axial vector strength, 
respectively. From these numbers we can conclude that we have reached the large- A limit 
already in the fp shell, supporting the conclusions of ||. We have also explained this 
saturation in qualitative terms. In 56 Ni and 100 Sn, where a major shell closure separates the 
spin-orbit partners, a larger effect is seen. This is caused by the first order contributions 
to the effective operator. The quenching stays nearly constant for energies up to some 60 
MeV in all cases. Therefore it is justified to speak about energy-independent quenching 
factors for beta decays in a given mass region. We also found that the second-order terms 
(in inverse nucleon mass) are relatively unimportant for most calculations. In particular, 
the uncertainties in the nuclear model calculations of the one-body transition densities mask 
these tiny corrections, generally a few percent at maximum. 

The quenching factors are used to extract the value of the ratio Cp/Ca at q = 100 MeV. 
In light systems 16 and 40 Ca our results indicate a small (of the order of few percent) 
quenching. In 100 Sn, we obtain an enhancement of the same order of magnitude. In 56 Ni a 
large quenching is seen. 

The next step in studies of effective operators is the inclusion of subnucleonic degrees of 
freedom in the evaluation of the different diagrams entering the definition of the effective 
operator. Especially we have in mind the A-isobars as an intermediate state. These states 
have essentially been neglected due the lack of a suitable AA interaction. We plan to 
extend our formalism to include such states through the use of a newly refitted nucleon- 
nucleon interaction which includes isobars as explicit degrees of freedom. This interaction 



[27] accounts for scattering data up to ~ 1 GeV in laboratory energy. The inclusion of 
meson exchange effects together with the effective transition operators is also a considerable 
task, not fully attacked yet. 

This work has been supported by the Academy of Finland under the Finnish Centre 
of Excellence Programme 2000-2005 (Project No. 44875, Nuclear and Condensed Matter 
Programme at JYFL). 
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TABLES 

TABLE I. Reduced nuclear matrix elements and the corresponding single-particle operators 
[without the lepton radial wave function G^{r)]. The j w {qr) are the spherical Bessel functions and 
^kwu are the vector spherical harmonics jj"6|| . The momentum operator for nucleons is p, and p v is 
the momentum of the neutrino. 



Matrix element 


Okwu 




[Owu] 


3w(qr)Y Q % u (r)5 wu 




[lwu] 


j w {qr)Y^ u {v,a) 




Ok wu ± 


[Owu±] 




j w {qr) ± aZ{m^/p u )j WTl {qr) 


Y(y£ u (r)8 W u 


[lwu±] 




j w {qr) ± aZ(m'^/p v )i WT i{qr) 


Y^ u (r,a) 


Okwup 


[Owup] 


ijw(qr)Y Q M u (r)a ■ P 5 WU 




[lwup] 


ijw(qr)Y^ u (r,p) 
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FIGURES 
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FIG. 1. Non- folded diagrams to second order in the interaction included in the evaluation of 
the effective operator. Hole states are represented by greek letters while particle states are given 

by roman letters. The operator itself is given by x in the various diagrams, while the wiggly 

lines are the nuclear G-matrix. Folded diagrams to second order in the interaction are included in 
the calculation but not shown here. 



23 




24 



0.85 



0.84 





1 1 1 
vector — i — 
axial — x— 
pseudoscalar 


i 


i 


1 1 










—¥c—~~~ 










1 






























- 
















i i i 


i 


1 


1 1 





10 20 30 40 50 60 70 80 90 100 

Q (MeV) 

FIG. 3. Renormalization of the vector, axial vector and pseudoscalar terms with 40 Ca as the 
closed-shell core, k = — 1 ( J = — > 1 transition, Air = no). 
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FIG. 4. Renormalization of the vector, axial vector and pseudoscalar terms with 56 Ni as the 
closed-shell core, k = — 1 ( J = — > 1 transition, Air = no). 
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FIG. 5. Renormalization of the vector, axial vector and pseudoscalar terms with 100 Sn as the 
closed-shell, k = — 1 ( J = — > 1 transition, A7r = no). 
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